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ABSTRACT 

Many Early-type galaxies have ionized gas emission in their centres that extends 
to scales of ^ lkpc. The majority of such objects are classified as low-ionization nu¬ 
clear emission regions (LINERs), but the nature of their ionizing source is still not 
clear. The kinematics associated with these gaseous structures usually shows devia¬ 
tions from a pure rotational motion due to non-gravitational effects (e.g. outflows) or 
to non-axisymmetric potentials (e.g. bars or tri-axial systems). This is the third of 
a series of papers that describes a sample of 10 nearby (d < 30 Mpc) and massive 
(cr > 200 km s -1 ) early-type galaxies observed with the Gemini Multi-Object Spec¬ 
trograph in Integral Field mode installed on the Gemini-South telescope. In paper 
II, we performed spectral synthesis to subtract the stellar components from the data 
cubes of the sample galaxies in order to study their nuclear spectra. In this work, we 
analyse the circumnuclear gas emission (scales of ~ 100 pc) of the sample galaxies 
and we compare the results with those obtained with PCA Tomography in paper I. 
We detected circumnuclear gas emission in seven galaxies of the sample, all of them 
classified as LINERs. Pure gaseous discs are found in three galaxies. In two objects, 
gaseous discs are probably present, but their kinematics are affected by non-Keplerian 
motions. In one galaxy (IC 5181), we detected a spiral structure of gas that may 
be caused either by a non-axisymmetric potential or by an outflow together with a 
gaseous disc. In NGC 3136, an ionization bicone is present in addition to five compact 
structures with LINER-like emission. In galaxies with a gaseous disc, we found that 
ionizing photons emitted by an AGN are not enough to explain the observed Ha flux 
along this structure. On the other hand, the Ha flux distribution and equivalent width 
along the direction perpendicular the gaseous disc suggest the presence of low-velocity 
ionized gas emission which seem to be related to the nuclear activity. We propose a 
scenario for LINER-like circumnuclear regions where a low-velocity ionization cone is 
formed by a collimating agent aligned with the gaseous disc. 

Key words: Techniques: imaging spectroscopy - galaxies: active - galaxies: elliptical 
and lenticular, cD - galaxies: ISM - Galaxies: kinematics and dynamics - galaxies: 
nuclei. 


1 INTRODUCTION 

The presence of ionized gas in the central regions of early- 
type galaxies (ETGs) has been studied by several works in 


the literature (|Phillips et al. 

1986 Kim 1989 Sarzi et al. 

2006 2010 

Yan & Blanton 

2012). However, the source 


of ionization is still not clear. Since most of this emis¬ 
sion is characterized as low-ionization nuclear emission re¬ 
gions (LINERs) ( Heckman||1980| Ho||2008 ), several mecha¬ 
nisms may be responsible for ionizing the gas: active galac- 
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tic nuclei (AGNs) - Ferland & Netzer 

< 

1983); 

Halpern & 

Steiner| (1983J; shocks - Heckman (1980 

i; 

Dopita & Suther- 


(pAGBs) - Binette et al 
Fernandes et al. (2011); 


(1994); Stasiriska et al. (2008); Cid 


Yan & Blanton (2012). The study 


of unresolved nuclear sources is very useful in detecting the 
presence of an AGN. Broad components in permitted lines 
and the presence of compact radio and X-ray cores are typ¬ 
ical signatures of nuclear activity in LINERs ( |Ho||2008 ). In 
addition, circumnuclear (scales of ~ 100 pc) and extended 
regions (scales of ~ 1 kpc) may also help to constrain the 
source of ionization in the central regions of ETGs. 
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Hubble Space Telescope (HST) observations of circum- 
nuclear regions of LINERs have shown that their kinematics 
may be more complex than pure rotational motion typical 
of gaseous discs ( Pogge et al.|2000 Walsh et al.|2008 ). More¬ 
over, only one LINER of the sample studied by |Pogge et al.| 
(2000) seems to have the same structure found in Seyfert 
galaxies, i.e. ionization cones and linear structures. The ex¬ 
ception is NGC 1052, which clearly shows an ionization cone 
structure aligned with a radio jet. With regard to photoion¬ 
ization, these authors suggest that a central UV source is 
responsible for the emission of the circumnuclear region, al¬ 
though they do not speculate on the nature of this source 
(it could be an AGN or hot stars). Walsh et al. (2008), using 
spectra obtained with the Space Telescope Imaging Spectro¬ 
graph (STIS), detected significant electron-density gradient 
in four LINERs and showed that, in most galaxies, the gas 
velocity dispersion peaks in the region dominated by the 
black hole and decreases in the outer regions. They argue 
that, in some galaxies, outflows may increase the gas veloc¬ 
ity dispersion and that the presence of virial motions close 
to the central supermassive black hole (SMBH) may either 
increase or decrease gas velocity dispersion in the nuclear 
region, depending on the emission line surface brightness or 
even on the inclination of the disc. Shields et al. (2007), 
comparing HST observations with data from the Palomar 
Survey (Filippenko h Sargent 1985} Ho et al.|199 7), showed 
that the gas emission is spatially distributed, with an exten¬ 
sion of ~ 100 pc, and that accretion-powered objects reveal 
a concentration of the emission in the nuclear region. 

In more extended regions, Sarzi et al. (2006) showed 
that the gas kinematics does not behave as a pure disc in 
a large fraction of the ETGs sample taken from SAURON 
integral field unit (IFU) data, although they proposed that 
the disturbances are probably caused by non-axisymmetric 
potentials (e.g. bars or tri-axial systems) on kpc-scales. In 
the same sample, Sarzi et al. (2010) noticed that the con¬ 
tribution of AGNs to the photoionization of LINERs is re¬ 
stricted to the nuclear region and that pAGB stars are likely 
responsible for the extended gas emission. The same findings 
were reported by Singh et al. (2013), using IFU data from 
the Calar Alto Legacy Integral Field Area (CALIFA) Sur¬ 
vey. Both Sarzi et al. (2010) and Singh et ah| ( |2013| ) highlight 
the fact that the association of LINER emission with nuclear 
activity may overestimate the population of AGNs in the 
Universe, since aperture effects may include regions domi¬ 


nated by extended photoionization sources (see e.g. Yan & 
Blanton 2012). In fact, |C id "F ernand es et al.| ( [2011| proposed 


that a comparison of the equivalent width (EW) of the Ha 
emission line with the [N II]/Ha ratio (the WHAN diagram) 
may be used to separate AGN from pAGB stars as the source 
of photoionization in LINER galaxies, even when aperture 
effects are present in the data (e.g. data from the Sloan Dig¬ 
ital Sky Survey - SDSS). Papaderos et al. (2013), also using 
data from the CALIFA survey, showed that LINER emis¬ 
sion from ETGs may be divided in two types: one where 
AGNs and pAGB stars are a dominating photoionization 
source in the nuclear and extended region, respectively, and 
the other where the leakage of Lyman continuum photons 
is important, leading to reduced emission-line luminosities 
and EWs. This could explain why galaxies with strong AGN 
emission in radio continuum are weak in the optical region. 
Furthermore, they argued that the relative distribution of 


gas when compared to the stars may be an important issue 
as well. Since in tri-axial systems the volume of the emitting 
gas is limited when compared to the stellar component, the 
nuclear EW is lowered in these cases, leading to an obser¬ 
vational bias against detection of AGN activity in gas-poor 
galaxy spheroids. Another possible scenario for weak optical 
emission and strong radio sources is that most of the energy 
released by the AGN may be mechanical and not radiative. 


This is the third of a series of papers that aim to char¬ 
acterize the nuclear and circumnuclear regions of a sample 
of 10 ETGs observed with the IFU installed on the Gemini 
South Telescope. In this paper, we intend to characterize 
the circumnuclear region of these 10 ETGs. In |Ricci et al.| 
(2014b| (hereafter paper I), we analysed this sample using 


Principal Component Analysis (PCA) Tomography (Heyer 
|V Schloerb|1997| Steiner et al.|2009| ). In |Ricci et al. (2014a) 
(hereafter paper II), we analysed the emission lines from 
the nuclear region after subtracting the stellar component 
from each spectrum of the data cubes. The results from both 
works showed that all 10 galaxies have typical LINER emis¬ 
sion and, in seven of them, we detected bona fide AGNs. We 
proposed that eight galaxies have circumnuclear emission 
and in seven of them the kinematics may be interpreted as a 
gaseous disc, while an ionization cone is likely to be present 
in one galaxy. However, using only the PCA Tomography 
technique, it is hard to detect kinematic perturbations in 
a gaseous disc by means of, let’s say, an outflow or even 
a non-axisymmetric potential. Thus, in this work, we will 
analyse the emission lines from the circumnuclear region af¬ 
ter subtracting the stellar component from the data cubes. 
It is worth mentioning that the CALIFA and the SAURON 
surveys also performed such a study, but in more extended 
regions (scales of kpc). However, both projects have poor 
spatial resolution to resolve the circumnuclear region. In 
the CALIFA project, the fibres of the IFU instrument have 
a diameter of 2.7 arcsec (Sanchez et al. 2012) and in the 
SAURON integral field spectrograph, the instrument allows 
a spatial sampling of 0.94 arcsec in the spatial low resolu¬ 
tion mode (Bacon et al. 2001). Although HST observations 
allow for such information in a scale from 10 to 20 pc, it 
is restricted to long-slit data. In the Gemini Multi-Object 
Spectrograph (GMOS) in IFU mode, the fibres have a diam¬ 
eter of 0.2 arcsec and the spatial resolution is limited by the 
seeing of the observations. Although limited to a field of view 
(FOV) of 3.5 arcsec x 5 arcsec, it allows one to characterize 
the circumnuclear region of the ETGs using 3D spectroscopy 
with superb spatial resolution for ground-based telescopes. 


Section [2] contains a brief description of the data. Sec¬ 
tion [3] describes the emission lines extracted from the cir¬ 
cumnuclear regions. In Section [4] we use BPT diagnostic 
diagrams to classify the emission type of the circumnuclear 
regions. In Section [5] we build maps of the radial velocity, 
velocity dispersion, Ha and [N II] fluxes and EWs, [N II]/Ha 
ratios and electron density for the whole FOV of the obser¬ 
vations. In Section [6] we discuss if the output from AGNs 
are enough to explain the circumnuclear emission of ionized 
gas. In Section [ 7 ] we discuss the presence of discs or ion¬ 
ization cones in the sample galaxies. In Section [8] we show 
images taken with the HST of some galaxies of the sample. 
Finally, in Section [9] we present a brief discussion of the re¬ 
sults and our main conclusions. In appendix [Aj we provide 
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comments for individual objects, summarizing the findings 
presented in papers I, II and III for each galaxy. 


2 BRIEF DESCRIPTION OF THE DATA 


The sample analysed in this paper was previously discussed 
in paper I. In short, it is composed by 10 massive ETGs 
(cr* > 200 km s- 1 ) with distances up to 30 Mpc. The galax¬ 
ies were observed with the Gemini South Telescope using 
the GMOS-IFU in one-slit mode (programmes GS-2008A- 
Q51 and GS-2008B-Q21). The data were reduced with the 
standard IRAF package for the Gemini telescopes. Bias, 
flat-fields and CuAr lamp images, in addition to spectropho- 
tometric standards, were acquired to properly correct and 
calibrate the data. The final product is a data cube with a 
FOV of 3.5 arcsec x 5 arcsec, where each spaxel has a size 
of 0.05 arcsec. Besides the basic reduction process, we ap¬ 
plied a Butterworth filter (low-pass) to each image of the 
data cubes with the aim of removing high frequency noise 
from the spatial dimension. A spectral low frequency instru¬ 
mental fingerprint was identified and duly removed with the 
PCA Tomography technique (Steiner et al. 2009, paper I). 
The data cubes were also corrected for the differential at¬ 
mospheric refraction (DAR) effect using an algorithm de¬ 
veloped by us. Finally, we deconvolved each image of the 
data cubes with the Richardson-Lucy deconvolution process 
( Ric hardsonj 1972~ Lucy|1974 ) with six iterations and assum¬ 
ing a Gaussian point spread function (PSF), whose FWHM 
values are equal to the seeing of the observations shown in 
paper I, where a more detailed description of the data cube 
reduction process and treatment of the noise is presented. 

We also performed a stellar population spectral synthe¬ 
sis in each spectrum of the data cubes. For this task, we 


used the starlight software (Cid Fernandes et al. 2005). 


The stellar population basis in our synthesis was described 


by Walcher et al. (2009). This basis uses the MILES library 
( Vazdekis et al.||2010 ) with solar abundances, whose fluxes 
were recalibrated pixel by pixel for non-solar abundances 
using the theoretical models of simple stellar populations 


(SSPs) of Coelho et al. (2007). The basis has 120 SSPs with 
a spectral resolution of 2.51 A ( Falcon-Barroso et al.||20lT 


Beih ori et al.||2011 ), ages between 3 and 12 Gyr with steps 
of 1 Gyr, abundances [Fe/H] = -0.5, -0.25, 0.0 and 0.2 and 
[a/Fe] = 0.0, 0.2 and 0.4. More details on the spectral syn¬ 
thesis performed on the data cubes of the sample galaxies, as 
well as the best-fitting stellar models for the spectra of the 
sample galaxies, are presented in paper II (Appendix A). Fi¬ 
nally, we subtracted the stellar features from each spectrum 
of the data cubes, resulting in gas cubes, i.e., data cubes with 
only the spectral lines related to the gaseous component. In 
paper II, we analysed only the spectra extracted from the 
nuclear region of the galaxies. In the present work, the entire 
FOV of the gas cubes is studied to extract information from 
the circumnuclear regions of the sample galaxies. 


3 CIRCUMNUCLEAR EMISSION LINE 
PROPERTIES FOR EACH SAMPLE 
GALAXIES 


We extracted representative spectra from the circumnu¬ 
clear regions of the sample galaxies. For each gas cube, we 
summed all spectra within a circular region. Since each gas 
cube has a distinct spatial dimension (due to DAR correc¬ 
tion) and also because of the position of the galactic nucleus 
of each object in the FOV, the extracted radius varies for 
each galaxy and it is limited by the smaller spatial dimension 
of the cube. These radii are presented in Table [2] The final 
circumnuclear spectrum of each galaxy is the result of the 
sum described above minus the respective nuclear spectrum 
that was extracted and analysed in paper II. 


The line profiles were fitted by means of the Gauss- 
Newton algorithm. It fits non-linear functions using the least 
squares method. More details about the fitting procedure 
may be found in paper II. The profiles were fitted with 
the sum of two Gaussian functions for each emission line. 
We fitted first the [S II]AA6716, 6731 doublet. For the other 
emission lines, the only free parameters of the fitting pro¬ 
cedure were the amplitudes of the Gaussian functions, once 
the gas velocity dispersion and the radial velocities of each 
Gaussian were constrained by the results found for the [S II] 
doublet (similar procedure was adopted by Ho et ak||199 7). 
Although different lines may have different profiles, since 
they arise from different regions, the kinematic constraint is 
very helpful to describe the emission lines in spectral regions 
where the signal-to-noise ratio is quite low. Even with these 
constraints, we did not detect the Hline of NGC 3136. The 
fitting procedure did indeed find a solution for H/3 in this 
galaxy, but the amplitude of the sum of the two Gaussian 
functions is ~ 1.3cr, where a is the standard deviation of the 
spectral region of this line. In ESO 208-G21, we did not de¬ 
tect the [O IJA6300 line to a level of 1.2cr. We did not detect 
any emission lines in the circumnuclear region of NGC 1399, 
NGC 1404 and NGC 2663. 

Figs [l] [2] and [ 3 ] show the fitting results for the [S II], 
Ha+[N II] and H/3+[0 III] lines, respectively. The Ha/H/3, 
[N II]/Ha, [S II]/Ha, [O I]/Ha and [O III]/H/3 line ratios, 
in addition to the flux calculated for the Ha line, are shown 
in Table [l] The errors presented in this table are purely sta¬ 
tistical and do not take into account systematic errors that 
may arise from an improper subtraction of the stellar com¬ 
ponent. However, it is worth mentioning that the absolute 
flux calibration is very uncertain for our data cubes, since 
the spectrophotometric standards and the galaxies were ob¬ 
served in different nights. Therefore, one should expect an 
error of about 30% or even higher in the absolute flux (or 
luminosity) measurements. The colour excess E(B-V), also 
shown in Table [l] was calculated with the observed Ha/H/3 
ratio using the extinction curve proposed by |Cardelli et al.| 
(1989) with Ry = 3.1 and assuming an intrinsic Ha/H/3 ra¬ 


tio for LINERs of 3.1 (Ferland V Netzer 1983 Halpern V 
Steiner 1 1983| . The Ha luminosities, corrected for reddening 


effects (except for NGC 3136, whose Ha luminosity is pre¬ 
sented without this correction), are shown in Table [ 2 ] For 
the electron densities n e of the circumnuclear regions, we 
used the [S II]A6716/[S IIJA6731 ratio. We calculated the 
values for n e with the NEBULAR package under the IRAF 
environment ( |Shaw V Dufou r 1995|, assuming T e = 10000K 
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and using the calculations proposed by De Robe rtis et al.| 
(1987). Due to the uncertainties associated with the [S II] 
lines ratio, we were able to assert only upper limits for the 
density of the circumnuclear regions of the sample galaxies. 
As a consequence, we estimated lower limits for the ionized 
gas mass of the circumnuclear regions, since this parameter 
is inversely proportional to n e (see Eq. 1 from paper II). 

The annular regions from which we extracted the spec¬ 
tra of the circumnuclear regions are spatially resolved. This 
allows us to calculate the filling factor of the gas within the 
volume of these cylindrical shells. The filling factor / is given 
by 


/ = 


L(Ha) 
e V n| ’ 


(i) 


where V is the volume of the cylindrical shell and e is the 
Ha line emissi vity. Assuming case B and T e = 10000 K, e 
= 3.84xl0 -25 (Osterbrock V Ferland 2006). For simplicity, 
we assumed that the height of the cylinder is equal to its ra¬ 
dius. Since we have only upper limits for n e , we calculated 
the lower limits for the filling factors of the circumnuclear re¬ 
gions of the sample galaxies. The results are shown in Table 

m 

The Ha luminosities of the circumnuclear regions are 
quite similar to those found for the nuclear regions (paper 
II, see also erratum). In five galaxies, the colour excess is 
consistent with zero, i.e., there are no significant reddening 
effects in the circumnuclear regions of these objects. Only 
for NGC 1380 a considerable E(B-V) was detected in this 
region. The upper limits of the circumnuclear densities are 
systematically lower than the nuclear densities. This result 
is in agreement with the proposals of Filippe nko &; Sargent | 
(1988), Ho et al. (1996) and Walsh et al.| ( |2008| b of a negative 
radial gradient for the ionized gas densities of the central 
regions of galaxies. For the ionized gas mass, the lower limits 
vary between 10 3 and 10 4 M©. 


4 DIAGNOSTIC DIAGRAMS 


Diagnostic diagrams were used to classify the circumnuclear 
emission of six galaxies of the sample (NGC 1399, NGC 
1404 and NGC 2663 have no circumnuclear emission and 
NGC 3136 has no detected H /3). We used the BPT dia¬ 
gram (Baldwin et al. i M. with the line ratios proposed by 
|Veilleux &; Oste rbr ock| ( |198 7). The values shown in Table [l] 
were inserted into the [N II]/Ha versus [O III]/H f3 and [S 
II]/Ha versus [O IIIJ/H/3 diagrams, together with the line 
ratios of the galaxies from the Palomar Survey ([Ho et al.| 
[1997] ). We also added the maximum theoretical limit for H II 
regions proposed by Kewley et al. (2001), the empirical sep¬ 


aration between AGNs and starburst galaxies proposed by 


Kauffmann et al. (2003) and the separation between LINERs 
and Seyferts proposed by Kewley et al. (2006). According to 


Kewley et al. (2006), the transition objects (TOs) are found 


in the region between the maximum theoretical limit for H 
II regions and the empirical separation between AGNs and 
starburst galaxies. The diagnostic diagrams are presented in 
Fig. El We found that the spectra of all six galaxies analysed 
with this tool have typical LINER line rations. 


5 MAPPING THE NUCLEAR AND 

CIRCUMNUCLEAR PROPERTIES OF THE 
NARROW LINE REGIONS 

One of the advantages of using data cubes is that we are 
able to obtain information of a given spectral region in two 
spatial dimensions. In this section, we present maps con¬ 
taining information about the ionized gas of the gas cubes 
of the sample galaxies, such as line fluxes, kinematics and 
electron densities. All maps were drawn by fitting Gaussian 
functions to the emission lines of each spectrum of a given 
gas cube. Before the fitting procedure, we subtracted the 
broad components of the Ha lines, related to the broad line 
regions (BLRs), from each spectrum of the gas cubes of the 
galaxies ESO 208 G-21, IC 1459, IC 5181, NGC 2663, NGC 
4546 and NGC 7097. We performed this subtraction because 
we are interested in analysing only the properties of narrow 
line regions (NLRs) of the galaxies, even on the position of 
the nucleus of each galaxy. To do so, we built data cubes 
whose spectra contained only the broad component of the 
Ha emission line. This was achieved by creating a normal¬ 
ized Gaussian image (i.e., the sum of all spaxels of this image 
equals one), where the FWHM is related to the PSF of the 
data cubes (papers I and II). Then, each spaxel was multi¬ 
plied by a spectrum containing only the broad component 
of the Ha line, which is also represented by a Gaussian func¬ 
tion. The radial velocity, flux and FWHM of the Gaussian 
related to the BLR were measured in the nuclear spectra of 
the galaxies and are shown in paper II. The result is a data 
cube that contains only the BLR of a given object, which 
has the same spatial and spectral dimensions as the respec¬ 
tive gas cube. The subtraction of the broad component from 
the spectra implies that the information taken from the Ha 
and from the [N II] emission lines are associated only with 
their narrow component (in fact, the [N II] doublet would be 
affected by the broad component of the Ha emission line), 
even in the nuclear region of these six galaxies. 

For each spectrum of the gas cubes, the fitting proce¬ 
dure is quite similar to that described in Section [3] Nev¬ 
ertheless, in regions of the FOV far from the nucleus, the 
signal-to-noise ratio is not enough for the procedure to fit 
more than one Gaussian function per emission line. Thus, 
each emission line was described by only one Gaussian func¬ 
tion. Moreover, the maps of NGC 1380 and NGC 2663 do 
not cover the whole spatial dimension of the gas cubes, since 
no emission lines exist in some regions in the FOV of both 
objects. In the case of NGC 1380, there is no gas emission 
in the upper region of the FOV. In NGC 2663, we were only 
able to map the emission lines of the region closest to the 
AGN. This is in accordance with the fact that we did not 
detect emission lines in the circumnuclear region of NGC 
2663. In NGC 1399 and NGC 1404, the emission lines are 
too weak and, therefore, no maps were constructed for either 
object. 

For the seven galaxies where we detected circumnuclear 
emission, and NGC 2663, we built maps of the gas kinemat¬ 
ics, flux and EW of the Ha and [N IIJA6583 emission lines. 
We also drew maps of the electron densities of the gas. In 
keeping with procedures in Section [3] we also used the [S 
II]AA6716, 6731 emission lines to calculate n e . Each map is 
discussed in the sub-sections below. 
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Figure 1. [S IIJAA6716, 6731 emission lines detected in the circumnuclear region of seven galaxies of the sample. The line profiles were 
fitted with a sum of two Gaussian functions, shown in blue and in green in the graphs above. For both [S II] lines, each Gaussian has 
the same kinematic parameters (radial velocity and velocity dispersion). The full fitted profile is shown in red. 
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Figure 2. Ha and [N IIJAA6548, 6583 emission lines detected in the circumnuclear region of seven galaxies of the sample. Both Gaussian 
functions were fitted with the same radial velocities and gas velocity dispersions as the ones fitted in the [S II] lines. 
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Figure 3. H (3 and [O IIIJAA4959, 5007 emission lines detected in the circumnuclear region of seven galaxies of the sample. Both 
Gaussian functions were fitted with the same radial velocities and gas velocity dispersions as the ones fitted in the [S II] lines. For a 
better visualization, we smoothed the observed spectra, although the fitting procedure was performed with the original data. 
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Table 1. Flux of the emission lines detected in the spectra of the circumnuclear regions of seven galaxies of the sample. The Ha flux is 
in units of 10“ 15 erg s -1 cm -2 . The colour excess E(B-V) parameter was calculated through the Ha/H/3 ratio. 


Galaxy name 

/(Ha) n 

(Ha/H (3) n 

E(B-V) 

[N II]/Ha 

[S II]/Ha 

[O I]/Ha 

[O III]/H/3 

ESO 208 G-21 

n±i 

2.68=1=0.60 

-0.14=1=0.22 

1.28=1=0.13 

1.26=1=0.14 


2.014=0.45 

IC 1459 

24±1 

3.09=1=0.19 

0.00=1=0.08 

1.80=1=0.08 

1.31=1=0.07 

0.17: 0.0 1 

1.704=0.13 

IC 5181 

13=1=1 

3.02=1=0.33 

-0.03=1=0.10 

1.10±0.10 

0.78=1=0.07 

0.07=1=0.01 

2.174=0.23 

NGC 1380 

6=1=1 

4.52=1=0.85 

0.36=1=0.18 

1.23=1=0.21 

0.77=1=0.14 

0.15=1=0.03 

1.274=0.25 

NGC 3136 

26=1=1 

- 

- 

1.62=1=0.10 

1.17=1=0.07 

0.19=1=0.01 


NGC 4546 

22=1=1 

2.58=1=0.23 

-0.18=1=0.09 

1.17=1=0.08 

0.95=1=0.06 

0 . 11 ± 0.01 

2.0.0, t 0.2 1 

NGC 7097 

5=1=1 

2.74=1=0.27 

-0.12=1=0.09 

0.98=1=0.09 

1.214=0.10 

0.144=0.02 

1.684=0.17 


Table 2. Physical characteristics of the circumnuclear component of the sample galaxies. Column (1): Ha luminosity in erg s — 1 . Column 
(2): Electron density in cm -3 . Column (3): Radius, in arcsec, of the outer region of the cylindrical shell used to extract the spectra from 
the circumnuclear regions of the galaxies. Column (4): Filling factor of the gas inside the cylindrical shell. Column (5): Ionized gas mass, 
in solar mass. In NGC 3136, the Ha luminosity is presented assuming E(B-V) = 0, since we did not detect the H (3 emission line in this 
object. 


Galaxy name 

log L(Ha) 

(1) 

n e 

(2) 

radius 

(3) 

f 

(4) 

Mlon 

(5) 

ESO 208 G-21 

38.584=0.29 

<153 

1.45 

>2.7xl0 —4 

> 1.8 X 10 3 

IC 1459 

39.324=0.14 

<90 

1.35 

>1.2x 10 —3 

>3.5 X 10 4 

IC 5181 

38.994=0.40 

<312 

1.5 

>5.8 X 10 —6 

>5.7 X 10 2 

NGC 1380 

38.754=0.24 

<658 

1.3 

>1.2x 10 —5 

>8.4 X 10 2 

NGC 3136 

39.254=0.27* 

<226 

1.25 

>1.9X 10 —4 

>6.7 X 10 3 

NGC 4546 

38.884=0.25 

<142 

1.45 

>3.4 X 10 —4 

>5.0 X 10 3 

NGC 7097 

38.764=0.14 

<100 

1.25 

>9.6 X 10 —4 

>8.7 X 10 3 


5.1 Gas kinematics maps 

For the kinematic properties, we fitted the [N II]AA6548, 
6583 and Ha emission lines of each spectrum of the gas 
cubes. The theoretical ratio [N II]A6583/[N II]A6548 = 3.06 
(Osterbrock & Ferland| [2006| was fixed for the fitting pro¬ 
cedure. We extracted the radial velocity and the velocity 
dispersion from each spectrum of the gas cubes of the sam¬ 
ple galaxies. All three Gaussians were assumed to have the 
same a and were separated in the Ha and [Nil] lines by 
their respective wavelength. In other words, there is only 
one value for the radial velocity and for the velocity disper¬ 
sion per spectrum. The velocity dispersion, which is related 
to the a of the Gaussian function and not to the FWHM, 
was corrected for the instrumental broadening effect as a 2 = 
a obs ~ a insti where the <Ji ns t values correspond to the spec¬ 
tral resolution of the data cubes (see paper I). The kinematic 
maps are shown in Figs. [5a| |5b| [5c| and |5d| for eight galaxies 
of the sample. 

In the radial velocity maps, we detected bipolar struc¬ 
tures along the circumnuclear regions in the eight analysed 
galaxies. The position angles (P.A.) of these kinematic bipo¬ 
lar structures were measured by fitting a linear function be¬ 
tween the points of maximum and minimum radial velocity 
values along this structure. The P.A. results are shown in 
Table [3] The errors of the P.A. were calculated by means 
of a Monte Carlo simulation, where we measured the P.A. a 
hundred times in the maps after adding Gaussian noise with 
cr equal to the error of the radial velocity in each spaxel. The 
uncertainties associated with the radial velocity are the sta¬ 
tistical errors taken from the fitting procedure of the emis¬ 
sion lines. 

In addition to the maps, we built ID profiles along the 
kinematic bipolar structure and along the perpendicular di- 


Table 3. Column (1): Galaxy name. Column (2): Radial veloc¬ 
ity, in km s —1 , measured on the nucleus of the galaxies after 
the redshift correction applied to the galaxies using the velocities 
presented in paper I. It is worth mentioning that these velocities 
were measured for the narrow component of the emission lines, 
since we removed the broad component of Ha from the gas cubes. 
Column (3): Position angles of the kinematic bipolar structures 
measured in the radial velocity maps shown in Figs. |5a| |5b| [5c] 
and |5d| Column (4) Position angles of the kinematic bipolar struc¬ 
tures measured in the maps of the blue and red wings of the Ha 
emission line shown in Fig. [ 7 ] 


Galaxy name 
(1) 

V nuc 
(2) 

P.A. Vr 

(3) 

P-A -Hoc 
(4) 

ESO 208 G-21 

13 

-66° ± 2 

-74° ± 5 

NGC 1380 

-50 

-169° 4= 2 

-169° ± 4 

IC 1459 

-65 

48° ± 1 

22° 4= 3 

NGC 7097 

-6 

-9° 4= 3 

-11° 4= 9 

IC 5181 

44 

-37° 4= 2 

-24° ± 4 

NGC 4546 

-12 

-130° 4= 2 

-150° ± 6 

NGC 2663 

66 

-110° ± 3 

-98° ± 6 

NGC 3136 

69 

78° 4= 2 

73° 4= 9 


rection of the same structure. In both cases, the profiles 
traverse the nucleus of the galaxies. The nuclear positions 
of the six galaxies with a BLR were estimated with the im¬ 
age of the red wing of the broad component of the Ha line. 
In both NGC 1380 and NGC 3136, we used the [O I]A6300 
emission line to define the nuclear position because these 
galaxies lack a BLR. For each galaxy, the ID profiles of 
the kinematic parameters are shown below their respective 
maps in Figs |5a| |5b| |5c| and |5d| Note that the velocity val¬ 
ues shown in Figs. |5a| |5b| |5c| and |5d| are presented relative 
to the velocity measured in the nucleus, whose values are 
shown in Table [3] It is worth mentioning that the velocity 
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[Nil]/Ha [S II] / Ha 


Figure 4. BPT Diagram with the line ratios measured in the circumnuclear spectra of six galaxies of the sample shown as filled magenta 
circles. The line ratios are presented in Table[l] According to the classification proposed by |Ho et al.| ( [1997| , the red crosses are Seyferts, 
green triangles are LINERs, hollow black circles are TOs and blue squares are H II regions. The thin brown line is the maximum starburst 
line proposed by Ke wley et al.| |2001| ), the dashed black line is the empirical division between H II regions and AGNs (Kauffmann et aT| 
2003) and the thick purple line is the LINER-Seyfert division suggested by |Kewley et al.| ( |2006| ). 


measured in the nucleus is relative to the redshifts of the 
sample galaxies shown in paper I. 

In paper I, we proposed that gaseous discs were present 
in six galaxies of the sample. In NGC 3136, using PCA To¬ 
mography, we detected the kinematics of an extended struc¬ 
ture, which is probably not associated with a gaseous disc. 
The radial velocity map of this object, shown in Fig. [5c| re¬ 
veals the overall kinematics of the structure. The P.A. of the 
kinematic bipolar structures detected with PCA Tomogra¬ 
phy (paper I) match the P.A. measured in the radial velocity 
maps of five objects. In NGC 3136, we found a difference of ~ 
20°, but the extended structure is quite irregular. In IC 5181, 
the difference between both measurements is ~ 16°. How¬ 
ever, a twist in the gas structure of this object was detected, 
which is discussed in Sections |5.4| and [7| In NGC 4546, the 
difference is ~ 19°. In this case, the ID profile along the 
kinematic bipolar structure suggests a non-Keplerian curve, 
with perturbations within a projected distance of 1 arcsec 
from the nucleus. Although IC 1459 has a difference of ~ 11° 
between both measurements of the P.A., the kinematic bipo¬ 
lar structure is affected by an outflow. This perturbation in 
the kinematics of IC 1459 has been previously mentioned by 


dispersion maps, we detected a peak 
in the nuclear region of all galaxies, except for NGC 3136, 
which has a velocity dispersion peak northward from the 
nucleus. In the circumnuclear regions of four galaxies, the 
variation of the velocity dispersion is the same in both par¬ 
allel and perpendicular directions of the kinematic bipolar 
structure. In the case of ESO 208 G-21, the velocity dis¬ 
persion falls more gently in the direction from the nucleus 
along the kinematic bipolar structure. In the case of NGC 
4546 and NGC 7097, the velocity dispersion radial gradient 
is shallower in the direction from the nucleus along the per¬ 
pendicular direction of the kinematic bipolar structure. In 
IC 5181, the velocity dispersion rises along the perpendic¬ 
ular direction of the kinematic bipolar structure. In NGC 
3136, a strong positive radial gradient is present northward 
from the nucleus. 

The gas kinematics suggest that pure gaseous discs are 


Cappellari et al. (2002 


In the gas velocity 


present only in ESO 208 G-21, NGC 1380 and NGC 7097. 
It is worth mentioning that in the first two objects, the 
gaseous disc is corotating with the stellar component (paper 
I). In NGC 7097, we detected a counterrotation between the 
gaseous disc and the stellar component (paper I). In NGC 
4546 and IC 1459, gaseous discs may be present but are af¬ 
fected by non-Keplerian effects, possibly outflows. In NGC 
3136, an ionization bicone is likely to be present. The cases 
of NGC 2663 and IC 5181 will be studied in Section 13 


5.2 Flux maps of the Ho and [N II]A6584 emission 
lines 

The same Gaussian functions used to extract kinematic in¬ 
formation from the gas structure of the sample galaxies were 
integrated in order to extract the flux of the Ho and [N 
II]A6583 emission lines. Figs. 


6 a 


and 6b show the results. 


The same figures present the ID profiles of the Ho flux maps 
along the parallel and the perpendicular directions of the 
kinematic bipolar structures. 

In the flux maps, we see a concentrated emission in the 


nuclear region of the galaxies. Indeed, Shields et al. (2007) 
showed that galaxies containing AGNs have a large concen¬ 
tration of the Ho flux in their nucleus. This is observed in 
the six galaxies containing a BLR and also in NGC 1380 and 
NGC 3136. This latter galaxy also has a slightly extended 
emission in the NE-SW direction. This feature was also re¬ 


ported by Goudfrooij et al. (1994). The extended emission 
is in the same position as the two compact objects detected 
with PCA Tomography applied to the spectral range 6250- 
6800 A (paper I). 

The ID profiles of the Ho flux reveal that, in six galax¬ 
ies of the sample, the circumnuclear regions are more intense 
along the kinematic bipolar structures. Moreover, in the di¬ 
rection perpendicular to the kinematic bipolar structures, we 
always detect emission and, from now on, we will call this 
“the low-velocity emission”, as can be clearly seen in Figs. 
[5a| [5b] [5c] and \Sd\ In NGC 3136, the ID profiles show that 
both directions have, approximately, the same behaviour. 
This is in agreement with the hypothesis that the geometry 
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Figure 5a. Left: Gas radial velocity map, in km s —1 , taken from the [N II]+Ha: emission lines. Right: Gas velocity dispersion map, in 
km s -1 . In the profiles, the hollow black circles were extracted in the direction parallel to the P.A. of the bipolar gas structures, while 
the filled green circles are related to the direction perpendicular to the P.A. of the bipolar structures. The red vertical lines delimit the 
FWHM of the PSFs of the data cubes. All the velocities shown above are relative to the nuclei of the galaxies. The radial velocities 
originally measured for the nuclear regions of the galaxies are shown in Table |3| 
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Figure 5b. Same as in Fig. |5a| 
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Figure 5c. Same as in Fig. |5a| 
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Figure 5d. Same as in Fig. |5a| 
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Figure 6a. Left: [N II]A6583 flux map, in erg s -1 cm -2 . Centre: Ha flux map, in erg s -1 cm -2 . Right: ID profile of the Ha flux 
maps. It is worth mentioning that the results shown above correspond to the NLR only, since we removed the broad components of Ha 
before the construction of the maps. The hollow black circles were extracted along the kinematic bipolar structures, while the filled green 
circles were extracted from the perpendicular direction. The red lines correspond to the PSF of the gas cubes. The blue dashed lines are 
the simple model of a spatially infinite, optically and geometrically thin gas structure with constant density and filling factor, which is 
photoionized by a central AGN. This model falls with R —2 . In each galaxy, the model was convolved with the PSF of the gas cubes. 
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Figure 6b. Same as in Fig. |6a| 
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of the circumnuclear region of this galaxy is not disc-like. In 
NGC 2663, the ID profiles in both directions are consistent 
with a point-like source. This is expected, since we did not 
detect any emission lines in the spectrum of the circumnu¬ 
clear region of this galaxy. 

In addition to the flux maps, we show in Fig. [7] images 
corresponding to the red and the blue wings of the narrow 
component of the Ha emission line (remembering that the 
Ha broad component was already subtracted). Again, we see 
bipolar structures in the circumnuclear region of the galax¬ 
ies. Moreover, the images shown in Fig.[7]are very similar to 
the tomograms related to the gas kinematics, which were ob¬ 
tained with PCA Tomography (paper I). We measured the 
P.A. of these structures and the results are shown in Table 
[3] The P.A. measured in these images are in agreement with 
the P.A. measured in their corresponding tomograms shown 
in paper I within 3a, except for IC 1459, whose difference be¬ 
tween both measurements is ~ 4a. The isophotes presented 
in Fig. [7] of ESO 208 G-21, NGC 1380 and NGC 7097 show 
narrow structures, which is compatible with the hypothesis 
of gaseous discs. In IC 1459, the outer isophotes are irregu¬ 
lar, probably affected by the outflows. This is also seen, at a 
lower intensity, in NGC 4546. Although the isophotes of IC 
5181 and NGC 2663 are disc-like, these cases are discussed 
in Section]?] For NGC 3136, the structure seen in Fig.[7]also 
resembles an ionization bicone. 


5.3 [N II] /Ha ratio maps 


The maps of the [N II]A6583/Ha ratio, shown in Figs. 8a and 


|8b| were drawn by dividing the flux maps of both emission 
lines. The ID profiles along the kinematic bipolar structures 
and along the low-velocity emission are also shown in Figs. 
I8alandl8bl 

The [N II]/Ha ratio has higher values in the nucleus of 
seven galaxies in the sample. In the circumnuclear region of 
galaxies ESO 208 G-21, IC 5181, NGC 4546 and NGC 7097, 


[N II]/Ha >1.0 along the kinematic bipolar structure. This 
is typical of LINER emission. In the low-velocity emission of 
these galaxies, [N II]/Ha < 1.0. In IC 1459, the ID profiles 
along both directions have similar behaviours, also seen in 


NGC 2663. 


In NGC 1380, the ID profile along the kinematic bipo¬ 
lar structure reveals that [N II]/Ha ~ 1.5 in a region located 
1 arcsec south from the nucleus. In the same direction, at 
1.8 arcsec from the nucleus, [N II]/Ha ~ 0.6, which is typ¬ 
ical of an H II region. Northward, also within a projected 
distance of 1.8 arcsec, [N II]/Ha ~ 0.8. Probably, a second 
H II region is to be found here. In NGC 3136, a few com¬ 
pact structures are seen in the map with [N II]/Ha > 1.8. 
One of the objects is likely to be the central AGN studied in 
papers I and II. In this galaxy, in the position of the second 
compact structure detected with PCA Tomography (paper 
I), we have [N II]/Ha ~ 1.6. 


5.4 Maps of the equivalent widths of the [N II] 
and Ha emission lines 

In order to construct the EW maps, each spectrum of the 
gas cubes was divided by the average of the continuum re¬ 
gions adjacent to the Ha absorption line in the spectrum 


taken from the stellar component. Subsequently, we inte¬ 
grated the [N II] and Ha emission lines in the same way de¬ 
scribed in Section [5^21 Th e EW(Ha) and EW([N II]) maps 
are presented in Figs. |9a| and |9b| 

In five of the eight galaxies shown in Figs. [9a| and [9b] 
the EW(Ha) is higher along the kinematic bipolar struc¬ 
ture than in the low-velocity emission. In NGC 2663, the 
ID profiles along both directions have the same behaviour. 
In all galaxies, the nucleus has the highest values of the EWs 
of both emission lines, except for NGC 1380, whose case is 
discussed below. In IC 5181, both EW maps reveal a spi¬ 
ral structure, similar to an integral sign. This structure is 
revisited in Section [7] 

In NGC 1380, the EW(Ha) map revealed an object lo¬ 
cated 1.8 arcsec south from the nucleus. In the same posi¬ 
tion, we have [N II]/Ha ~ 0.6. The EW([N II]) map does 
not show any structure in this region. These results reinforce 
the hypothesis that an H II region is present at this posi¬ 
tion. At 1.8 arcsec northward from the nucleus, we detected 
an extended object in the EW(Ha) map. The [N II]/Ha ra¬ 
tio values in this position also suggest an H II region. In 
addition to these star-forming regions, the EW([N II]) map 
shows a compact object located at 0.7 arcsec south from the 
nucleus, in the same position where [N II]/Ha ~ 1.5. In the 
EW(Ha) map, this object is not seen. This may be the con¬ 
sequence of a second LINER located in the central region of 
NGC 1380. This discussion is revisited in Section [8] and also 
in a future paper. 

Maps of the EW of both lines of NGC 3136 show at least 
three extended structures. One corresponds to the structure 
detected in the flux maps of both lines, which may be re¬ 
lated to two compact objects detected with PCA Tomog¬ 
raphy (paper I). The other two structures, one located 1.4 
arcsec southeastward from the nucleus and the other one 
located 1.5 arcsec southwestward from the nucleus have [N 
II]/Ha ~ 1.5 and ~ 1.3, respectively. Both ratios are typical 
of LINERs. 


5.5 Density maps 

Following procedures in Section [ 3 ] we used the [S II] A6716 
/ [S II] A6731 ratio to calculate the electron density in each 
spaxel of the gas cube. The [S II] fluxes were extracted using 
the same procedure described in Section [5.2| However, only 
the amplitudes of the Gaussian functions were assumed as 
free parameters. Radial velocity and velocity dispersion were 
constrained to the values found for the [N II] and Ha emis¬ 
sion lines. The electron density was calculated by means of 
the relation proposed by Proxauf et al. (2014J, assuming T e 
m 10000K. Note that in Section [3] we used the NEBULAR 
package under the IRAF environment to calculate n e . As we 
used this package to estimate the nuclear densities in paper 
II, we were able to compare the densities of both the nuclear 
and circumnuclear regions. On the other hand, [Proxauf et ah] 
(2014) proposed an empirical numerical function to calcu¬ 


late n e , which was easier to implement in the code that we 
used to extract the flux from the emission lines. Within the 
errors, both calculation methods provide the same results. 
The density maps are shown in Figs. |10a] and |10b] 

In the outer regions of the FOV, the signal-to-noise ratio 
of the [S II] lines is quite low. This may be the cause of sev¬ 
eral artefacts seen in Figs. |10a] and |10b] since the measure- 
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Figure 7. RGB images corresponding to the blue and red wings of the Ha emission lines. It is worth mentioning that the results shown 
above correspond to the NLR only, since we removed the broad components of Ha before the construction of the maps. The image of 
the blue (red) wing minus the image of the red (blue) wing is shown in blue (red). The tip of the PSF corresponding to the flux maps 
in Figs. |6a| and |6b| is shown in green. We also added isophotes for each structure, which are shown as white contours. 


ment of the density is very uncertain. Thus, it is quite hard 
to interpret these maps. Still, some things may be worth 
mentioning. For instance, in NGC 2663, the density seems 
to be constant along the FOV analysed in this paper. In 
NGC 3136, the densest regions do not coincide with any 
object mentioned before, not even the nucleus. In ESO 208 
G-21, two structures are seen in the respective density map, 
wherein none of them are in the same position as the AGN. 


6 IS AGN PHOTOIONIZATION ENOUGH TO 
EXPLAIN THE CIRCUMNUCLEAR 
EMISSION OF THE SAMPLE GALAXIES? 


We want to check if the circumnuclear emission of the sam¬ 
ple galaxies can be explained by photoionization from their 
AGN. As a first approach, we will adopt the same model 
as Sarzi et al. (2010) and Singh et al. (2013), that is the 


case where the flux distribution of the Ha; emission line is 
caused by a point-like ionizing source embedded in an in¬ 
finitesimally thin gas disc of constant density and filling fac¬ 
tor being optically thin to Lyman continuum radiation. In 
this simplified model, the Ha emission line flux is expected 
to fall as R -2 , where R = 0 sets the position of the AGN. In 
fact, the flux distribution of Ha follows the strength of the 
radiation field emitted by an AGN. It is a simplistic model, 


since an infinitesimally geometrically thin and optically thin 
gas distribution may not apply for real ETGs. Moreover, the 
assumptions of a constant density and filling factor do not 
hold for real objects. For instance, a negative 
distribution was observed for some galaxies 

). Given the 

the model, it is quite helpful to verify if the photoionization 
caused by an AGN is enough to explain the circumnuclear 
emission. For each galaxy, this model was convolved by the 
PSF of the respective gas cube. The results are shown to¬ 
gether with the ID profiles of the observed Ha emission line 
flux in Figs. [6a| and 1 6b | 

Along the low-velocity emission, the observed profiles 
are somehow consistent with the model results. In some 
cases, as ESO 208 G-21, IC 1459 and NGC 4546, the ob¬ 
served ID profile along this direction is steeper than the 
model predictions. One possible reason is that these objects 
have a negative radial gradient of densities. Another hypoth¬ 
esis is that the nebula is not optically thin around the re¬ 
gions closest to the AGN. In NGC 1380 and NGC 3136, the 
observed ID profile along the low-velocity emission is shal¬ 
lower than the model predictions. Probably, the presence 
of other compact ionization sources apart from the central 
AGN of both NGC 1380 and NGC 3136 may be affecting 
the observed ID profiles in both directions (see Section |8|. 
It is worth mentioning that in NGC 3136, the kinematic 
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Figure 8a. Left: [N II]/Ha ratio maps. Right: ID profiles of the [N II]/H a ratio. It is worth mentioning that the results shown above 
correspond to the NLR only, since we removed the broad components of Ha before the construction of the maps. The hollow black circles 
correspond to the observed ratio along the kinematic bipolar structure, while the filled green circles are related to the observed ratio 
along the low-velocity emission. The vertical red lines delimit the FWHM of the PSFs of the data cubes. 
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Figure 8b. Same as in Fig. |8a| 
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Figure 9a. Left: [N II] EW maps, in A. Centre: Ha EW, in A. Right: ID profiles of the Ha EW. The hollow black circles were extracted 
along the kinematic bipolar structures, while the filled green circles were extracted from the low-velocity emission. The vertical red lines 
delimit the FWHM of the PSFs of the data cubes. 
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Figure 9b. Same as in Fig. |9a| 
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Figure 10a. Left: Density maps in cm -3 . Right: ID profiles of the density along the parallel and perpendicular directions of bipolar gas 
structure. The hollow black circles were extracted along the kinematic bipolar structures, while the filled green circles were extracted 
from the low-velocity emission. The vertical red lines delimit the FWHM of the PSFs of the data cubes. 
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bipolar structure is probably related to an ionization bicone 
while for the other objects we found a disc-like geometry. 
In NGC 7097, the results from the model match the ob¬ 
served ID profile along the low-velocity emission. However, 
it does not mean that this region of NGC 7097 follows the 
very simplistic assumptions of the model. For this galaxy, it 
is possible that the gas distribution along the low-velocity 
emission is not geometrically thin, which would result in a 
shallower profile caused by projection effects. Indeed, this 
geometry may be true for the other galaxies of the sample, 
except for NGC 3136 whose geometry is probably associ¬ 
ated with an ionization bicone. On the other hand, the ID 
profiles of the flux distribution of Ha along the kinematic 
bipolar structures are much shallower that the model pre¬ 
dictions. Although projection effects may also be happening 
here, this is probably not enough to account for the missing 
photons emerging from the gas structure (see e.g., the re¬ 
sults of Sin gh et al.| 2013). Therefore, other photoionization 
sources may be important to explain part of the circumnu- 
clear emission of these galaxies. 

Another way to test the importance of AGNs in pho- 
toionizing the circumnuclear regions is by means of pho- 
toionization models. To do so, we used the nebular spectral- 
simulation program CLOumQ We assumed a plane-parallel 
geometry, a power-law continuum f u oc v a , with a ~ -1.5, 
which is typical for LINER-like AGNs ( Ho|2008] ), a lower cut 
in the energy of the ionizing photons of 27 eV (photons with 
less energy than 27 eV have already been absorbed by the 
most internal regions of the nucleus), a filling factor of 10 -3 
and a gas density of 200 cm -3 . We also used the ionization 
parameter U for the calculations, which is the ratio between 
the ionizing photon density and the gas density. Since an 
AGN is a compact source, we may calculate U as 


( 2 ) 


v = rh = QW 

n e Aiir 2 n e c^ 

where n 7 is the ionizing photon density, c is the speed of 
light and Q(H) is the number of ionizing photons (i.e. with 
energy > 13.6 eV) emitted by the AGN. It is previously 
known that LINER-like AGNs have log U ~ -3.5 (Ferland 
Netze"r||1983| Halpern &; Steiner||1983| |Ho||2008 ). For 


calculations, we found that -3.5 < log U < -3.4 produces 1.81 
< [O III]/H/5 < 2.98, [N II]/Ha - 1.10, [S II]/Ha - 0.77, 
and [O I]/Ha ~ 0.11. These results match the observed line 
ratios shown in Table[l] For log U — -3.4 and n e — 200 cm -3 , 
an AGN must have n 7 = 0.079 cm -3 in order to produce a 
LINER-like spectrum. 

In order to verify how far an AGN may ionize the cir¬ 
cumnuclear region, we have to estimate Q(H) for the sample 
galaxies. With the nuclear Ha luminosity (paper II and erra¬ 
tum) and assuming that the nebulae are radiation bounded, 
it is possible to estimate the number of ionizing photons 
emitted by the AGNs. This is calculated as 


Q(H) = 


L Ha a B (H°,T) 


(3) 


Ill'll a OtHct(H° ,T) ’ 

where as is the recombination coefficient for all levels of 
the H atom, a Ha is the recombination coefficient for the Ha 
line, h is Planck’s constant and VHa is the wave frequency 


1 Calculations were performed with version 07.02.00 of CLOUDY, 


last described by Ferland et al. (1998). 


Table 4. Number of ionizing photons emitted by the AGN of the 
sample galaxies. 


Galaxy name log(Q(if)) 


ESO 208 G-21 

51.15 

fl¬ 

0.22 

IC 1459 

51.74 

it: 

0.13 

IC 5181 

51.06 

± 

0.39 

NGC 1380 

50.61 

± 

0.24 

NGC 2663 

51.88 

± 

0.44 

NGC 3136 

50.66 

± 

0.19 

NGC 4546 

51.14 

± 

0.40 

NGC 7097 

51.38 

± 

0.13 


Table 5. Column (1): galaxy name. Column (2): Projected dis¬ 
tances at which the ionizing photons emitted by the nucleus of 
the sample galaxies have a density of 0.079 cm -3 . This value for 
the ionizing photon density is produced by an AGN with log U ~ 
-3.4, n e ~ 200 cm -3 that, according to photoionization models, 
is able to explain the line ratios shown in Table [5] Column (3): 
same projected distances presented in column (2), but in pc. 


Name 

r 

r 


(arcsec) 

(P c ) 

(1) 

(2) 

(3) 

ESO 208 G-21 

0.85 

70 

IC 1459 

1.06 

139 

IC 5181 

0.53 

64 

NGC 1380 

0.43 

38 

NGC 2663 

1.25 

163 

NGC 3136 

0.34 

40 

NGC 4546 

0.84 

69 

NGC 7097 

0.61 

92 


of the photon related to the transition of the Ha line. We 
used as — 2.59x 10~ 13 cm 3 s -1 and ana — 9.39xl0 -14 cm 3 
s _1 (Osterbrock &; Ferland|2006 ). Table [Z] shows the results 
of the observed Q(H) of the sample galaxies. Assuming log 
U — -3.4 and n e — 200 cm -3 , we estimated the projected 
distance where the observed Q(H) for the sample galaxies 
produce the value for n 7 predicted by the photoionization 
model. The results are shown in Table [5] Note that all pro¬ 
jected distances are smaller than 1.25 arcsec. Figs. [6a| and 
[6b] suggest that all sample galaxies have ionized gas beyond 
1.25 arcsec along the gaseous discs. Also, if we assume that 
an AGN is photoionizing the region with a projected dis¬ 
tance of 2 arcsec, then to maintain n 7 = 0.079 cm -3 the 
number of ionizing photons should be 10 times higher than 
the observed values for the sample galaxies. These results 
reinforce the need for alternative photoionization sources. 


7 DISCS OR OUTFLOWS? THE CASES OF IC 
5181 AND NGC 2663. 

In three galaxies of the sample, we propose that the kine¬ 
matic bipolar structures are pure gaseous discs. In NGC 
1380 and ESO 208 G-21, the P.A. of the gaseous disc and 
the P.A. of the stellar disc are the same (paper I) and, prob¬ 
ably, the gaseous disc had an internal origin. Although the 
gas and the stellar component are counterrotating in NGC 
7097 (paper I, C aldwell et al.|1986 ), the radial velocity map 
and ID profiles (Fig. s in addition to the image of the red 
and blue wings of the Ha emission line, suggest that the gas 
component has a disc-like geometry. In IC 1459 and NGC 
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4546, a gaseous disc may exist, but they are affected by non- 
Keplerian motions, probably outflows (see Section 5.1). In 
NGC 3136, an ionization bicone seems very likely. 

An interesting case is IC 5181. In this galaxy, the planes 
defined by the kinematics of the stellar and gas components 
are perpendicular (paper I, Pizzella et al.|[20l3 ). Although 
an external gas accretion is able to produce a gaseous disc 
perpendicular to a stellar disc (Sarzi et al. 2006), which is 
the scenario proposed by Pizzella et al. (2013), an outflow 
is likely to be present in this object. The velocity dispersion 
map shows that lower values are seen along the kinematic 
bipolar structure, which is compatible with a geometric thin 
disc. On the other hand, the maps of EW(Ha) and EW([N 
II]) show a spiral structure, similar to an integral sign, in 
the inner regions (r < 0.7 arcsec). Two possible scenarios 
are suggested. One is that outflows are happening in regions 
very close to the AGN. The other possibility is that the gas 
is moving through a non-axisymmetric potential. This would 
create a concentration of gas in certain regions that would 


be outlined by the EW images of the emission lines (Sarzi 
et al.|2006 ). 


Another curious case is NGC 2663. The radial veloc¬ 
ity map reveals a bipolar structure (Fig. [5c]). The isophotes 
shown in the map of the red and blue wings of the Ha emis¬ 
sion line (Fig. [7| are regular. Although both properties are 
typical of a gaseous disc, these kinematics are seen within 
the unresolved region of the nucleus. This may be the reason 
why this galaxy has high values for the FWHM of the Gaus¬ 
sian profiles representing the emission lines of the spectrum 
extracted from the nucleus (paper II). Although we may not 
discount the possibility that NGC 2663 has a gaseous disc, 
the hypothesis of an ionization bicone is also plausible here. 
This is the only case for which the AGN has enough photons 
to photoionize the region along the bipolar structure (see 
Fig. 6b). The spatial extent is smaller when compared to 
the other galaxies (the galaxies IC 1459, IC 5181 and NGC 
7097 have, approximately, the same spatial scale as NGC 
2663). Although these two characteristics could be applied 
to a gaseous disc, an ionization bicone seems very probable 
in this galaxy. 


8 HST OBSERVATIONS OF THE SAMPLE 
GALAXIES 

Observations made with the HST may provide additional 
information to the results obtained with the gas cubes, es¬ 
pecially because of the superb spatial resolution of the HST 
images. ESO 208 G-21, NGC 1380, IC 1459 and NGC 3136 
have public image^] available in the HST data archive. Im¬ 
ages taken with the V and I filters and with a narrow band 
filter on the Ha+[N II]AA6548,6583 emission lines were taken 
from the data archive (see Table [6] for more information). 
Both Ha/I and Ha/V ratios highlight the ionized gas emis¬ 
sion of the galaxies. We are mostly interested in analysing 
compact regions of line emission as well as the central core 
of the bulges. In order to emphasize this, we performed a 

2 NGC 1380: GO 10240, PI - Aaron Barth; IC 1459: GO 6537, 
PI - Tim de Zeeuw, GO 5454, PI - Marijn Franx; NGC 3136: GO 
6822, PI - Paul Goudfrooij; ESO 208 G-21: GTO/ACS 9293, PI 
- Holland Ford 


Table 6. Summary of the HST images shown in Fig. m 


Galaxy name 

Instrument 

Filters 

ESO 208 G-21 

ACS/WFC1 

F658N(Ha + [N II]), F814W(/-band) 

IC 1459 

WFPC2/PC 

F555W( V-band), FR680P15(Ha + [N II]), 
F814W(/-band) 

NGC 1380 

ACS/HRC 

F555W( V-band), F658N(Ha + [N II]), 
F814W(/-band) 

NGC 3136 

WFPC2/PC 

F547M( V-band), FR680P15(Ha + [N II]), 
F814W(/-band) 


wavelet decomposition of the HST images. We used the “a 


trous” algorithm for the wavelet decomposition (Starck & 
Murtagh 2002). In Fig. [TT] we show only the first two wavelet 


components, W0 and Wl, which are related to the high spa¬ 
tial frequencies of the images. The colour V - I maps the 
intrinsic extinction of the galaxies and, eventually, points 
out the existence of a non-stellar continuum, as seen in IC 
1459. In Fig. ED we show HST images of the four galaxies 
mentioned above, with the same FOV and the same spa¬ 
tial orientation as the gas cubes. The comparison between 
the HST images and the data cubes analysed throughout 
this paper and also in Papers I and II took as reference the 
nuclear region of the galactic bulges. 

In IC 1459, the ionized gas emission image reveals an 
extended structure, probably associated with the gaseous 
disc, as shown in the Ha/V image (Fig.[TT]). The high spatial 
frequency wavelet component of the ionized gas emission 
shows a compact object in the central region of the galaxy, in 
the same position that a blue compact object is seen in the V 
-1 colour image. This is probably related with the non-stellar 
continuum of the AGN of IC 1459. This is in agreement with 
the detection of the featureless continuum of the AGN in the 
eigenspectrum that was associated with the AGN in paper I. 
A spiral structure, seen in the V - I colour image, probably 
formed by dust, is present along the southwest-northeast 
direction. Many of these HST results were already discussed 
by Carollo et al. (1997) and Cappellari et al. (2002). 

In NGC 1380, the ionized gas emission image shows 
three compact objects along the same FOV of the gas cubes. 
Two objects are separated by 0.2 arcsec along the southwest- 
northeast direction and are not spatially resolved in the gas 
cubes. However, the EW([N II]) map, shown in Fig. [9a] re¬ 
veals two objects separated by 0.7 arcsec. The separation 
between the two closest objects and the third object in the 
HST image is ~ 0.7 arcsec. In the region of all three ob¬ 
jects, [N II]/Ha > 1.2. These three objects could probably 
be classified as LINERs. A further analysis of the nature of 
these three objects will be discussed in a future paper. The 
H II regions, proposed in the sections above, are not seen in 
the HST images. With the V - I colour image, we detected a 
dust disc which extends along the entire north-south direc¬ 
tion of the FOV, to the west of the nucleus, which may be 
related to the external disc of NGC 1380, since this galaxy is 
classified as SO (see paper I). A smaller dust structure, with 
the same orientation of the external disc, is located in the 
nuclear region. This is in agreement with the high E(B-V) 
values found for the nucleus (paper II). 

In ESO 208 G-21, an intense ionized gas emission is 
seen in the central region of the galactic bulge. We also 
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NGC 1380 

Ha/I I-W0 + W1 V-l 




Figure 11. HST data archive of NGC 1380, NGC 3136, ESO 208 G-21 and IC 1459. In all galaxies, the images have the same FOV and 
spatial orientation as the gas cubes. For some images, we show the sum of the wavelet components WO and Wl. This components are 
associated with the high spatial frequencies of the images. The magenta circles mark the nuclear region of the bulges of the galaxies. In 
V - I colour images, the dark regions are associated with regions with large extinctions. 


detected an extinction along the east-west direction in the 
FOV, slightly to the south of the nucleus. 

In NGC 3136, the ionized gas emission is complex. We 
confirmed that a compact emission arises from the centre of 
the bulge. A weaker and slightly extended emission is seen 
northeast of the centre. Both objects are in the same posi¬ 
tion as the compact objects detected with PCA Tomography 
(paper I). We also noticed that the extended structures de¬ 
tected southeast of the nucleus in both the EW(Ha) and 
the EW([N II]) maps are related to two compact objects in 
the HST Ha/I image. The structure located southwest of 
the nucleus is at the border of the FOV of the gas cube, 


but it is extended in the Ha/I HST image (see also Fig. Al 


in Appendix A8). Moreover, the entire FOV related to this 
HST image seems to be filled by a diffuse gas emission. The 
V - I colour image shows an extended structure with a high 
extinction located southeast of the nucleus which is in the 
same position as the structure in redshift of the ionization 
bicone. This result indicates that this structure in redshift is 
the far side of the bicone, since it is probably located behind 
a dust screen that causes the extinction detected in the V - 
I colour image. 


9 DISCUSSION AND CONCLUSIONS 

In this paper, we analysed the circumnuclear emission of 
eight from a sample of 10 galaxies. The other two galaxies, 
NGC 1399 and NGC 1404, did not have circumnuclear emis¬ 
sion detected with PCA Tomography (paper I). After per¬ 
forming stellar spectral synthesis and subtracting the stellar 
components, an attempt to search for some type of ionized 
gas emission was successful only for the nuclear region of 


both objects (paper II). As mentioned in Section [3j we used 
PCA tomography (paper I) to detect circumnuclear emis¬ 
sion of ionized gas as kinematic bipolar structures, where 
one pole is redshifted and the other is blueshifted relative to 
the centre of the galaxies. The tomograms associated with 
these emissions resemble the maps corresponding to the blue 
and red wings of the Ha emission line, shown in Fig. [T| The 
kinematic bipolar structure is also seen in the radial veloc¬ 
ity maps of the galaxies. When compared to the nucleus, 
the circumnuclear region is less dense, but the lower limits 
for the ionized gas mass are always higher than the ion¬ 
ized gas mass in the unresolved nuclear region. In paper I, 
we proposed that six galaxies of the sample have gaseous 
discs. The results presented in this paper revealed that only 
ESO 208 G-21, NGC 1380 and NGC 7097 seem to have pure 
gaseous discs, i.e., their motions are apparently Keplerian. 
Furthermore, in NGC 1380 and ESO 208 G-21 we detected 
a corotation between the gaseous disc and the stellar com¬ 
ponent, as opposed to NGC 7097, whose gas structure and 
stellar component are counterrotating. However, this only 
means that, while the gaseous discs of ESO 208 G-21 and 
NGC 1380 may have an internal origin, the gaseous disc of 
NGC 7097 probably has an external origin, since the stel¬ 
lar radial velocity of this galaxy is slower than the gas ra¬ 
dial velocity (Caldwell et al. 1986 Ricci 2013). Moreover, 
the maps of the red and the blue wings of the Ha emis¬ 
sion line of these three objects have narrow and well-defined 
isophotes, in agreement with the geometry of a gaseous disc. 
On the other hand, it is very likely that IC 1459 and NGC 
4546 also have gaseous discs, but their radial velocities are 
affected by non-Keplerian motions. This is probably associ¬ 
ated with outflows within the central region of these objects. 
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For IC 1459, the non-Keplerian effects along the P.A. of the 
main gas structure were also detected by |Cappellari et al.| 


with an exponent -a yields a projected power-law intensity 


( 2002 ). In NGC 4546, Sarzi et al. (2006) observed twists 


along the radial velocity map and also a spiral structure in 
the ionized gas distribution, similar to an integral sign. It 


is worth mentioning that both Sarzi et al. (2006) and Gal- 


letta (1987) suggested a counterrotation between the stellar 


and the gas components for this galaxy, while we proposed 
in paper I that both structures are not in the same plane. 
However, the twists were detected in a region beyond the 
spatial coverage of the data cubes analysed in this paper. 
The same spiral structure was detected in IC 5181, but in a 
spatial scale covered by the FOV of the GMOS data cube. 


Sarzi et al. (2006) suggested that this spiral structure de¬ 


rives from non-axisymetric potentials in the central region 
of some ETGs. Nevertheless, we do not discard the hypoth¬ 
esis of an outflow in the most internal regions (<0.7 arcsec) 
of IC 5181. Even for NGC 4546, the scenario of an outflow 
for the observed twists may be true. In the case of NGC 
3136, the kinematic results together with the V - I colour 
map from the HST observations suggest the presence of a 
bicone, where the structure in redshift is the far side of the 
bicone, since an extinction was detected in the V - I colour 
image of the same region, which may indicate that this side 
of the bicone is behind a dust screen. In NGC 2663, we pro¬ 
pose two hypotheses: a bicone or the presence of two AGNs. 

Emission line ratios taken from the spectra of the cir- 
cumnuclear regions suggest LINER-like emission. Although 
the nuclear regions have higher [N II]/Ha ratio on the maps 
of seven galaxies, the kinematic bipolar structures also have 
typical line ratios of LINERs along their entire spatial ex¬ 
tension detected in the gas cubes. Only in NGC 1380 we did 
detect line ratios typical of H II regions in two positions of 
the gaseous disc. 

LINER-like spectra may be produced by several ioniza¬ 
tion sources (AGNs - |Ferland fc Netzer 1983| Halpern & 


Steiner|1983| shocks -|He ckman|1 980| pAGB stars -|Binette 


et al.|1994| |Stasinska et al.|2008| |Cid Fernandes et al.|201l[ 


Yan & Bl anton|2012 ). BLRs, which are a typical signature of 
an AGN, were detected in six galaxies of the sample. How¬ 
ever, the results presented in this paper suggest that a single 
AGN is not enough to photoionize the gaseous discs. Fol¬ 


lowing the work done by Sarzi et al. (2010) and Singh et al. 


(2013), we used a simple model to predict the Ha emission 


line flux distribution in an infinitisemaly thin and optically 
thin gas structure, with constant density and filling factor, 
which is ionized by an AGN. Since the simplistic assump¬ 
tions of this model do not hold, caution should be exercised 
when comparing the predictions from it with our measure¬ 
ments along the gaseous disc: The gas is not optically thin 
to Lyman continuum radiation, the Ha flux distribution is 
therefore not expected to decrease as R ~ 2 . The same situa¬ 
tion applies when the gas density and filling factor decreases 
with radius. On the other hand, the flux distribution will be 
flatter than R~ 2 under the assumption of a geometrically 
thick disc, since one has to take into account the integrated 


light over the line-of-sight. However, Singh et al. (2013) have 
shown that profiles that result from the projection of an ax- 
isymmetric disc with a radially declining gas density still do 
not explain the extended emission from their sample galax¬ 
ies. In an isotropically emitting-spheric symmetric volume, 
a luminosity density distribution decreasing as a power law 


profil e with an exponent -a+1 (see, e.g., Papaderos et al. 
2013). Therefore, a R~ 2 dilution of the radiation from a 


central point source (the AGN) will yield an intensity pro¬ 
file of the form R~ a with a — 1.0. We estimate a = 1.1 to 
1.7 along the gaseous discs of our sample galaxies. In this 
scenario, the AGN may be playing an important role in pho- 
toionizing the gaseous discs, however one would have to take 
into account the negative radial gradient of the gas density 
and also dust inside the gaseous disc. It seems very likely 
that an alternative ionization source has to be taken into 
account in addition to the AGNs. 

Photoionization models also show that the AGNs do 
not emit enough photons to explain the emission from the 
gaseous discs, since the line ratios of these structures are 
described by an ionization parameter of log U ~ -3.4. Table 

[5] shows that the sample galaxies have ionized radii beyond 
the limit that an AGN with such a value for U is able to 
produce. However, if we assume that the AGNs of the sample 
galaxies are able to photoionize regions with radii which are 
twice the values of the ionized radii calculated for log U = 
-3.4, then the ionization parameter in these farthest radii 
should be log U = -4.0. The same model we used in Section 

[ 6 ] but assuming log U = -4.0 produces a [N II]/Ha ratio of 
~ 0.74. Since the low-velocity emission of the galaxies has 
[N II]/Ha < 1.0, it seems reasonable that the AGN may be 
photoionizing these regions. However, it is worth mentioning 
that these photoionization models are also subject to errors. 
For instance, for the NLR radius calculation, we used the 
number of ionizing photons, which is proportional to the 
Ha luminosity. However, the uncertainty associated with the 
absolute flux calibration of the data cubes is quite high, 
certainly ^ 30%. Also the gas density and filling factor are 
very uncertain for the circumnuclear regions. Thus, caution 
is needed when analysing these results. 

Another result that supports the scenario where pho¬ 
toionization by a central AGN does not hold for the gaseous 
disc but seems to be important along the low-velocity emis¬ 
sion is presented in Fig.[l2] On the left of this figure, we show 
the relation between the typical EW(Ha) for the gaseous 
disc and the EW(Ha) of the nuclear region, while on the 
right we show the relation between the typical EW(Ha) for 
the low-velocity emission and the EW(Ha) of the nuclear re¬ 
gion. We notice a correlation between the nuclear EW(Ha) 
and the EW(Ha) of the low-velocity emission, while no 
correlation is seen between the nuclear EW(Ha) and the 
EW(Ha) in the gaseous disc. Given these results, we pro¬ 
pose a scenario where the low-velocity emission is related 
to ionization cones which arise in LINER-like AGNs, but 
the ionizing photons emitted by the AGN are collimated by 
some agent that is aligned with the gaseous disc. A scheme 


of this model was presented by Ricci et al. (2011) for NGC 
7097 (see their Fig. 5). 

One hypothesis for the collimation agent is that a nu¬ 


clear toms, within the context of the unified model (An- 
tonucci|199 3), exists in these galaxies and it is aligned with 


the gaseous discs. If this is true then, assuming that the inner 
radius of the toms is given by Rdust = O. 4 L 45 5 pc (Nenkova 
et al. 2008), where L 45 is the bolometric luminosity in units 
of 10 4t> erg s-\ we found that the Rdust is smaller than 
1 pc for all galaxies of the sample with a BLR (the bolo¬ 
metric luminosities of the sample galaxies are shown in the 
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Figure 12. Left: EW(Ha) along the gaseous discs versus nuclear EW(Ha). Right: EW(Ha) along the low-velocity emission (that is, 
extended emission perpendicular to the kinematic bipolar structure) versus nuclear EW(Ha). A correlation is seen only along the low- 
velocity emission. This may be an indication that AGN photoionization is important only along this region, while other ionization sources 
may be present along the gaseous discs. 


erratum of paper II). However, low-luminosity AGNs may 


Elitzur & Ho 


2009). According to the model proposed by NGC 1380 using X-ray data analysed by Su & Irwin (2013) 


lack a nuclear torus (Elitzur & Shlosman 2006 Ho 2008 


We calculated the Ha luminosity that would be emitted by a 
gas photoionized by the hot ISM in the galaxies IC 1459 and 


Elitzur &; Ho 2009), the BLR and the torus are formed by 


the same mechanism, i.e. winds emerging from the accretion 
disc. The boundary between the existence and the lack of a 
BLR/Torus structure occurs when an accretion disc around 
a SMBH has a very low radiative efficiency, which is the case 
for LINER-like AGNs ( |Hop 008). For a BLR/Torus to exist, 
L45 > 6.1 x 10 ~ 17 (REdd )~ 2 1 where REdd is the Eddington 
ratio. All sample galaxies are above this limit by a factor 
of ~ 100, except for NGC 1399 and NGC 1404. Thus, a 
nuclear torus is probably present in the eight galaxies anal¬ 
ysed in this work. If the torus is aligned with the observed 
gaseous discs, then the nuclear ionizing radiation would not 
reach this structure, but it would be able to photoionize the 
regions that are perpendicular to the gaseous disc, i.e. the 
low-velocity emission. On the other hand, the assumption 
that the nuclear torus is always aligned with the gaseous 
disc may not hold. However, if the orientation of the nuclear 
torus is randomly distributed when compared to the gaseous 
disc orientation, then we should expect that the gas distri¬ 
bution along the the low-velocity emission is not ionized for 
some galaxies. However, we can draw no firm conclusions on 
this issue because our sample is statistically incomplete and 
small. Alternative hypotheses for the collimation agent in 
the sample galaxies are that the inner radius of the gaseous 
discs is optically thick or some structure somehow related 
to the gaseous discs is acting as a shield against the ionizing 
photons from the AGN. 

Below, we propose some hypotheses for the ionizing 
sources: 

(i) AGN : In paper II, we confirmed that an AGN is 
present in, at least, seven galaxies of the sample, since in 
six of them we detected a BLR and NGC 1399 is a radio 
galaxy (Shurkin et al. 2008). In NGC 1380, NGC 1404 and 
NGC 3136, an AGN is likely to exist (papers I and II). 

(ii) Hot interstellar gas : ETGs have a hot interstellar 
medium (ISM), with a temperature of ~ 10 7 K. The ISM 
emits basically only X-ray photons and may be a source of 
photoionization for the circumnuclear region of the galaxies. 


We found that the ratio between the calculated and the ob¬ 
served Ha emission line luminosities is 0.02 for IC 1459 and 


0.10 for NGC 1380. Kim (1989), using a similar procedure, 


found similar results for three other ETGs with LINER-like 
nuclei. Therefore, the combined contribution of a hot ISM 
and an AGN does not seem to be enough to explain the Ha 
emission from the gaseous discs. 

(iii) Shocks : in NGC 3136, the presence of several struc¬ 
tures whose spectra have line ratios typical of LINERs with 
a bicone increases the probability of shocks in the central 
region of this galaxy. Supernovae remnants (SNR) or su¬ 
perwinds may produce a mechanical source of energy which 
may cause shocks along this galaxy. 

(iv) HII region contamination : In NGC 1380, we detected 
H II regions at a projected distance beyond ~ 1 arcsec from 
the nucleus. Although starbursts occupy a specific region in 
BPT diagrams, we should not discount the hypothesis that 
the light from the AGN is contaminating the H II region 
emission, which would characterize the observed spectrum 
as a transition object. 

(v) Hot low-mass evolved stars (HOLMES): |Binette et al. 
(1994) proposed that part of the ionizing photons observed 


in LINERs could be provided by pAGB stars. Several works 
in the literature ( Ho||2008| Stasiriska et al.||2008 |Eracleous 


et al.|20i0 Sarzi et al. |2Q1Q Cid Fernandes et al. 2011) have 

supported this idea. The values found for the EW(Ha) of the 
gaseous discs are smaller than 3A for all sample galaxies (see 
Fig. 12). According to the WHAN diagram (Cid Fernandes 
et al. || 2QTT ), this is an indication of photoionization caused 
by HOLMES. However, EW measurements may be affected 
by dilution effects caused by the stellar component in triaxial 
ETGs ( Papaderos et al.|2013 ). As a result, the observed EW 
would mimic HOLMES photoionization even in cases where 
photoionization caused by an AGN dominates the emission, 
which may be the case of ESO 208 G-21 (nuclear EW(Ha) 
< 3A, see Fig. Sal. Nevertheless, one of the consequences 
of the HOLMES scenario is that the EW(Ha) should be 
constant along the extended structure that is photoionized 
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by an extended source (Sarzi et al.| 2010). In the ID proflle 
of EW(Ha) along the gaseous discs, shown in Figs. 9a and 
|9bl we see this trend in ESO 208 G-21, IC 1459, IC 5181, 
NGC 4546 and NGC 1380. In NGC 7097, the EW(Ha) val¬ 


ues decrease until a projected distance of ~ 1.5 arcsec and 
then rise again. Thus, it is very likely that these old hot 
stellar populations are contributing significantly along the 
kinematic bipolar structures. However, it is not clear why 
they contribute only along the kinematic bipolar structures 
and not in the perpendicular direction, since the stellar discs 
are not always coupled to the gas motions and the bulge of 
the galaxies is almost isotropic. In order to verify the im¬ 
portance of HOLMES, we used CLOUDY to estimate the 
ionizing photon density emitted by HOLMES. Assuming a 
blackbody with T=150000K, n e = 200 cm -3 , a filling fac¬ 
tor of 10 -3 , a chemical composition of 2.5 times the solar 
metallicity and a log U ~ -3.3, we found [N II]/Ha ~ 0.91, 
[S II]/Ha - 0.80, [O I]/Ha - 0.08 and [O III]/H f3 ~ 3.36. 
Only the [O III]/H [3 ratio does not match the observations. 
This is probably related to the fact that we assumed a black- 
body radiation instead of using a spectrum of a real pAGB 
population. 


It is worth mentioning that several works in the liter¬ 
ature have discussed the importance of HOLMES in very 
extended regions, i.e., in scales of kpc (Sarzi et al. 2010 


Cid Fernandes et al.|2011||Yan &; Blanton|2012||Singh et al.| 

2013| . Given the superb spatial resolution of the GMOS data 
cubes, we are able to propose that it is very likely that pAGB 
stars are also very important in spatial scales of ~ 100 pc. 
Furthermore, it is important to emphasize that the greatest 
contributions of HOLMES is along the disc-like structure. 
This discovery was only made possible by the use of data 
cubes. 

Below, we summarize the main findings of this work: 


• We confirm the presence of kinematic bipolar structures 
in the circumnuclear region of the eight galaxies of the sam¬ 
ple where PCA Tomography revealed this feature (paper I). 
When comparing the radial velocity maps to the tomograms 
shown in paper I, the direction of the gas kinematics as well 
as the P.A. of the structures are in agreement within the 
errors. 

• In three galaxies, we detected pure gaseous discs in 
their circumnuclear regions. In other two objects, gaseous 
discs are probably present, but are probably affected by non- 
Keplerian motions (e.g. outflows). 

• In IC 5181, we detected a spiral structure of gas in the 
very central region of the galaxy. This may be the result of 
an outflow combined with a gaseous disc or even that the gas 
kinematics is the consequence of non-axisymetric potentials. 

• In six galaxies, we detected low-velocity emission, per¬ 
pendicular to the gaseous discs. 

• We detected an ionization bicone along the circumnu¬ 
clear region of NGC 3136. We also detected five compact 
objects in this galaxy. They all have LINER-like emission. 
Although we did not ascertain the nature of the photoion¬ 
ization source of these objects, it is possible that they are 
associated with shocks (SNR or superwinds). 

• We propose that the circumnuclear emission along the 
disc-like structures are not photoionized by the central AGN. 
The best candidate for the photoionization source is hot, old 
stellar population (e.g. pAGB stars). 


• Along the low-velocity emission, the emission seems to 
be consistent with photoionization by AGNs. 

• We propose a scheme for LINER-like circumnuclear re¬ 
gions where a low-velocity ionization cone, perpendicular to 
the disc, is formed by the collimation of the ionizing photons 
by a collimating agent aligned with the gaseous disc. 
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APPENDIX A: COMMENTS ON INDIVIDUAL 
OBJECTS 

Here, we summarize the main results from this work and 
also from papers I and II for each galaxy of the sample. 
Only for NGC 3136 a more detailed analysis is presented 
below. It is worth mentioning that the scenario proposed for 
the emission of ionized gas from the circumnuclear region 
discussed in Section [9] is common for all sample galaxies, 
except for NGC 1399, NGC 1404 and NGC 2663. 


A1 ESO 208 G-21 

Both nuclear and circumnuclear spectra have emission lines 
typical of LINERs. An AGN is responsible for the nuclear 
emission, since a broad component is seen in the Ha line 
(paper II). The circumnuclear gas component is probably a 
pure gaseous disc, as suggested by the radial velocity map 
(Section [5.1 1 and also by the fact that this gaseous structure 
is corotating with the stellar component (paper I). We also 
detected a reddening of the stellar component in the nuclear 
region caused by the interstellar medium (ISM) of this object 
(paper I). A peculiarity of ESO 208 G-21 is the fact that the 
density map reveals two structures that are not related to 
the AGN. However, a deep study of the nature of these two 
structures is beyond the scope of this paper. 


A2 IC 1459 


IC 1459 has a well-known LINER-like AGN (see references 
cited in paper I), thus this object was used to validate the 
results obtained with PC A Tomography in paper I. We con¬ 
firmed the presence of the AGN, since a broad component 
in the Ha line is present in the nuclear spectrum (paper II). 
Furthermore, we detected a featureless continuum from the 
AGN with PCA Tomography (paper I). The gaseous com¬ 
ponent in the circumnuclear region also has a LINER-like 
spectrum. If it is a gaseous disc, then the kinematics are 
affected by non-Keplerian motions (e.g. outflows, see Sec¬ 
tion 5.1) and it is counterrotating with respect to the stellar 


component (paper I). Cappellari et al. (2002), with long-slit 
data, reached the same conclusions for the kinematics of the 
circumnuclear region. 


A3 IC 5181 


This object also has LINER-like spectra for both nuclear 
(paper II) and circumnuclear regions. A broad component is 
seen in the Ha line (paper II). The gaseous component in the 
circumnuclear region has a spiral form, similar to an integral 
sign. This result may be related to an outflow in the nuclear 
region or to the fact that the gas is moving through a non- 
axisymmetric potential, which causes a gas concentration 
in certain regions (see Section [Tf . It is worth mentioning 
that the kinematics associated with this gaseous structure 
are perpendicular to the kinematics related to the stellar 
component (paper I). The same finding was also reported by 
Pizzella et al. ( |2013| . A reddening of the stellar component 
in the nuclear region, caused by the ISM of IC 5181, was 
also detected (paper I). 


A4 NGC 1380 


Both the nuclear and circumnuclear regions have LINER- 
like spectra. Although we did not detect a broad component 
in any permitted lines from the nuclear spectrum, the HST 
image shown in Fig. [Tl] reveals three compact sources in the 
central region of this galaxy; one of the objects is located 
at the centre of the bulge and the others are separated by 
0.2 and 0.7 arcsec from this central source. The EW([N II]) 
map in Fig. [9a] shows an extended structure that may be 
related to these compact objects, which indicates that they 
probably have LINER-like spectra. A deep analysis of these 
three structures will be performed in a future work, although 
we believe that the central object is probably related to an 
AGN. We also detected two H II regions near the centre 
of NGC 1380, seen in the EW(Ha) map (Fig. [9a] see the 
discussion in Section 5.4). A pure gaseous disc is present 


in NGC 1380 and it corotates with the stellar component 
(paper I). The reddening of the stellar component in the 
nucleus, caused by the ISM of this object, was also detected 
(paper I). 


A5 NGC 1399 


This galaxy has a weak emission of ionized gas in the nuclear 
region (paper II). O’Connell et al. (2005) detected evidence 
of an AGN in a far-UV spectrum taken with HST. Shurkin 
et al. ( |2008| showed that jets originated in an AGN created 
a cavity in the hot gas from the central region of NGC 1399. 
Thus, the emission lines detected by us in paper II are prob¬ 
ably related to this AGN. In paper II, we showed that most 
of the energy released by the AGN is mechanical. No circum¬ 
nuclear emission of ionized gas was detected. Also, we did 
not detect any kinematics related to the stellar component. 


A6 NGC 1404 


This is the only object where an AGN was not previously 


mentioned in the literature, although Grier et al. (2011) de¬ 


tected a point-like x-ray source with Lo^skeV ~ 10 u erg 
s _1 . They mentioned that this measurement is probable as¬ 
sociated with a low-luminosity AGN. We detected the Ha 
and [N II] emission lines in the nuclear spectrum. The bolo- 
metric luminosity estimated for this AGN, using the Ha line, 
matches the X-ray luminosity (paper II). We propose that 
NGC 1404 has a LINER-2 AGN. No circumnuclear emission 
of ionized gas is seen in this object. We detected a rotation 
of the stellar component with PCA Tomography (paper I). 
The reddening caused by the ISM was detected in a V - I 
colour image from the HST and also in eigenvector 2 of PCA 
Tomography applied to the 4325-6800 A spectral range. 


A7 NGC 2663 

The gas emission detected in NGC 2663 is restricted to the 
nuclear region. It has a LINER-like spectrum with a broad 
component in the Ha line. Although the emission is unre¬ 
solved, a kinematic signature of the gaseous component is 
detected, but it is probably related to an outflow rather than 
a disc. No additional ionization source is needed to explain 
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Table Al. Ha luminosities in erg s —1 and the gas density in 
cm -3 of the five objects detected in the central region of NGC 
3136. The objects marked with a * did not have their Ha lumi¬ 
nosities corrected for reddening effects. 


Object 

Spatial coordinates 

log L(Ha) 

n e 

1 (nuc) 

(0.0,-0.2) 

38.83=1=0.19 

267 

2 

(-1.0,0.0) 

38.41=1=0.08* 

130 

3 

(-0.3,-1.4) 

38.604=0.24 

< 100 

4 

(0.2,-1.1) 

38.344=0.08* 

120 

5 

(1.3,0.7) 

37.924=0.34 

140 


the gas emission of NGC 2663. When it comes to the stel¬ 
lar structure, no kinematic signature and also no reddening 
caused by an ISM were detected. 


A8 NGC 3136 


Along this paper and also in paper I, we have proposed that 
NGC 3136 may contain other photoionization sources, in 
addition to the central AGN. PC A Tomography applied to 
the spectral region between 6250 and 6800 A of the data 
cube detected two compact objects (hereafter objects 1 and 
2, see Fig. |A1| ). One of them is redshifted and the other 
is blueshifted relative to the centre of NGC 3136. More¬ 
over, both EW(Ha) and EW([N II]) maps revealed three 
extended structures, wherein the position of one of these 
structures corresponds to the objects 1 and 2. The struc¬ 
ture southeastward from the nucleus is related to objects 3 
and 4, while the structure southwestward from the nucleus 
was associated with object 5. In this section, we will try to 
ascertain the nature of these objects 

We extracted spectra of all five objects assuming a 
Gaussian PSF with a FWHM equal to the seeing of the 
observation (paper I) and with an amplitude given by the 
flux in the spaxel where the respective object is centred. The 
coordinates relative to the central region of the FOV of each 
object is shown in Table [Al] In objects 2 and 4, the spectral 
regions corresponding to both [O IJAA6300, 6363 and H/3 
emission lines are dominated by noise. In object 3, we did 
not detect the [O IJAA6300, 6363 doublet. Where Kf3 is de¬ 
tected, the Ha emission line luminosity shown in Table [A~T| 
was corrected for reddening effects as described in Section 
|3j The emission lines fluxes are shown in Table lA2l 

Objects 1, 3 and 5 were plotted in BPT diagrams [N 
II]/Ha versus [O III]/H (3 and [S II]/Ha versus [O III]/H£, 
which are shown in Fig. |A2| They may be classified as LIN- 
ERs. Objects 2 and 4 did not have H/3 detected, but their [S 
II]/Ha and [N II]/Ha are typical of LINERs. None of the ob¬ 
jects have a broad Ha component. However, the kinematic 
bipolar structure detected for this galaxy is probably asso¬ 
ciated with a bicone emerging from an AGN. Indeed, the V 


- I colour image taken with HST (see Fig. 11) shows extinc¬ 
tion in the region of the bicone that is redshifted relative to 
the centre of NGC 3136. This may indicate that the struc¬ 
ture in redshift is the far side of the bicone, behind a dust 
screen. Object 1 is in the centre of NGC 3136 and is proba¬ 
bly associated with the AGN of the galaxy. According to the 
PC A Tomography results (paper I), the Ha + [N II] lines of 
object 2 are redshifted when compared to the same lines of 
object 1. One possible explanation is that both objects are 


LINER-like AGNs and they are rotating around their centre 
of mass. Another possibility is that objects 2, 3, 4 and 5 are 
associated with shocks. The presence of the bicone acting 
as the source of mechanical energy and the fact that shocks 
may produce LINER-like spectra ( Heckman|1 98Q; Dopita V 


Sutherland 2003; Ho 2008) make this hypothesis quite likely. 


H II regions may also play a role in the central region of NGC 
3136. However, their [N II]/Ha and [S II]/Ha ratios are too 
high to be produced by young stellar populations, so one of 
the sources mentioned above (AGN or shocks) would have 
to coexist with this starburst regions in order to explain the 
high emission from low ionization elements. 


A9 NGC 4546 


We detected emission of gas in both the nuclear and cir- 
cumnuclear regions. The nuclear spectrum is in the border 
between LINERs and Seyferts and it also has a broad com¬ 
ponent in the Ha line. The circumnuclear gas emission has a 
LINER-like spectrum. The structure of this gaseous compo¬ 
nent is probably related to a gaseous disc, but non-Keplerian 
motions are seen in the radial velocity map (Section |5.1| . We 
also detected a rotation of the stellar component which is 
perpendicular to the motion of the gaseous structure. |Sarzi| 
et al. (2006) found that both structures are counterrotating. 


However, they also reported that the gaseous component 
has a spiral form. Although this twist is located in a re¬ 
gion beyond the FOV of our observations, our results are in 


accordance with those obtained by Sarzi et al. (2006). 


A10 NGC 7097 

This galaxy also has LINER-like emission in both the nu¬ 
clear and circumnuclear regions. In the nuclear spectrum, 
a broad component in the Ha line is present. For the cir¬ 
cumnuclear region, the gaseous structure is probably a pure 
gaseous disc, although it is counterrotating with respect 
to the stellar component. A reddening of the nuclear stars 
caused by the ISM was also detected. It is worth mentioning 
that the scenario for the circumnuclear gas emission pro¬ 
posed in this work (Section |9| w as previously suggested for 
this object in Ri cci et ah] (2011). 
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Figure Al. Left: Combined JHK infrared images of NGC 3136 obtained with 2MASS jSkrutskie et a l. |2006| . Right: HST image of the 
(Ha+[N II])/I ratio. This is the same image shown in Fig. but with a FOV of 7 x 7 arcsec. The FOV of the data cube is delineated 
in green. The numbers correspond to the objects detected in the central region of NGC 3136. Object 1 is the central AGN proposed in 
papers I and II. 


Table A2. Flux of the Ha emission line in 10 15 erg s 1 cm 2 


and line rations for the five objects detected in NGC 3136. 


Object 

/( Ha)„ 

(Ha/H/3)n 

E(B-V) 

[N II] /Ha 

[S II]/Ha 

[O I]/Ha 

[O III]/H/3 

1 (nuc) 

8.1±0.3 

3.39=1=0.57 

0.08=1=0.16 

1.82=1=0.09 

0.97=1=0.05 

0.19=1=0.03 

1.66=1=0.30 

2 

3.7±0.2 



1.70=1=0.15 

1.32=1=0.12 



3 

2.7±0.2 

4.27=1=0.92 

0.31=1=0.21 

1.57=1=0.17 

1.16=1=0.12 


2.12=1=0.51 

4 

3.2=1=0.2 



1.49F0.15 

1.03=1=0.09 



5 

1.2=1=0.2 

3.21=1=1.01 

0.03F0.31 

1.40F0.32 

0.91=1=0.17 

0.25=1=0.06 

1.56F0.47 


i-1—.—... 1 —>—.. 



0.01 0.1 1 0.01 0.1 1 
[Nil]/Ha [S ll]/Ha 


Figure A2. BPT Diagrams with the results of the five objects found in the central region of NGC 3136. The line ratios are presented in 
Table |A2| According to the classification proposed by |Ho et al.| (1997), the red crosses are Seyferts, green triangles are LINERs, hollow 
black circles are TOs and blue squares are H II regions. The thin brown line is the maximum starburst line proposed by |Kewley et ah] 
( |2001| ), the dashed black line is the empirical division between H II regions and AGNs (Kauff mann et ahj2 003 ) and the thick purple line 
is the LINER-Seyfert division suggested by |Kewley et al.| ( |2006| >. All five objects are in the LINER region of the diagram. 






















